2, 2008 6:20 WSPC/INSTRUCTION FILE OpenProb 



OPEN PROBLEMS IN HYDRODYNAMICAL APPROACH 
TO RELATIVISTIC HEAVY ION COLLISIONS 



T. Kodama, T. Koide, G. S. Denicol and Ph. Mota* 

We discuss some open problems in hydrodynamical approach to the relativistic heavy 
ion collisions. In particular, we propose a new, very simple alternative approach to the 
relativistic dissipative hydrodynamics of Israel and Stewart. 
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1. Introduction 

The phase transitions in strongly interacting bulk matter predicted by the quan- 
tum chromodynamics (QCD) should manifest their existence in different physical 
scenarios such as the evolution of inhomogeneities in the early universe, structure 
of compact stars, spectra of particles from nuclear collisions at ultra-relativistic 
energies, etc. In particular, the relativistic heavy ion collisions are the unique pos- 
sibility of observing such phase transitions in laboratories, permitting us to extract 
the bulk properties of the strongly interacting matter at extremely high tempera- 
ture and energy density. 

As we see from many talks in this conference, the over-all "picture" of the new 
states of the matter at extreme condition achieved in the laboratories (mainly from 
SPS and RHIC) is now being configured after more than two decades when the first 
project on relativistic heavy ion physics started^. Yet the nature of the transition 
from the hadronic phase to the QGP phase is still to be clarified quantitatively. Of 
course, in addition to the experimental data, advances in theoretical studies such 
as lattice QCD calculations also enriched understandings of the properties of the 
strong interacting matter^''^. 

Among many signals of QGP, the one which brought a new insight is the emer- 
gence of collective flow in the final state of exploding particles. In the non-central 
collisions, an asymmetric distribution of energy density is created in the first in- 
stant of the collision. In the hydrodynamical image, the driving force to expand 
the system is proportional to the pressure gradient, so that this asymmetry with 
respect to the direction of the largest pressure gradient should increase as a func- 
tion of the transverse momenta of particles. Quantitatively, such anisotropy can be 
expressed in terms of the coefficients of Fourier series of the azimuthal distribution 
of particles. We define the elliptic flow for a given transverse momentum window 
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as 




angles of the particle and of the impact parameter vector with respect to a some 
space-fixed coordinate system. In practice, the determination of these coefficients 
from the experimental data is not trivial since the reaction plane is not given a 
priori. In the hydrodynamical calculations, of course the reaction plane is given 
from the beginning. 

The finite positive value of this coefficient V2 is referred to as the elliptic flow 
and, from the point of view of hydrodynamics, it is sensitive to the initial pressure 
gradient of the system. The above expected behavior of V2 was clearly observed 
in RHIC data, and the values of observed values approach to those calculated 
from the ideal hydrodynamics. At the same time, for lower (SPS) energies, the 
experimental values are lower than the ideal fluid values (See Figs.l and 2). This 
fact was interpreted as the emergence of the new state of the matter which flows 
almost as an ideal fluid, while in the hadronic phase the matter suffers from the 
collisional viscosity**. 



Fig. 1. Elliptic flow at SPS energies as function of transverse momentum, compared to the simple 
ideal fluid hydrodynamical calculation. The observed values of elliptic flow are below the hydro 
values. Figure adapted from the presentation of M. Gyulassy, see this volume 

In addition to the behavior of the elliptic flow phenomena, hydrodynamical 
description is found to be very successful in many aspects of SPS and RHIC data^, 
thus establishing the general physical scenario of strongly interacting matter in 
the process of relativistic heavy ion collisions. On the other hand, this remarkable 
success of the hydrodynamics opens several important questions and understanding 
of these questions may lead us to some new physical concepts. In this paper, we 
discuss some of these questions. 




^fr, . . ■ ,1 .... ...1 . ., .. . .. ,.1 ., ■ . .. I... . . ... t ....... I. 




2, 2008 6:20 WSPC/INSTRUCTION FILE OpenProb 



Open Problems in Hydrodynamical Apprach to Relativistic Heavy Ion Collisions ) 3 



Au+Au; = 200 GeV; Mid-raptdity 




Transverse Momentum (GeV/c) 

Fig. 2. Elliptic flow at RHIC energies as function of transverse momentum, compared to the simple 
ideal fluid hydrodynamical calculation. Figure adapted from the presentation of M. Gyulassy, see 
this volume 

2. Relativistic Hydrodynamics 

The hydrodynamical description in nuclear and particle physics is not new^. How- 
ever, it is also true that to establish a solid theoretical foundation for the ap- 
plication of hydrodynamical description to microscopic systems has always been 
difficult. This is basically because the starting assumption of the validity of local 
thermodynamical equilibrium in hydrodynamics is not at all trivial to be justified 
for microscopic systems. 

For the sake of later discussion, let us here summarize briefly the basic structure 
of the relativistic hydrodynamics*. Let T^'' be the energy- momentum tensor of the 
matter. Then the dynamics of the system should obey the conservation law of this 
tensor, 

S^T"'' = 0. (2) 

For the case of a perfect fluid in local equilibrium, we can write 

T^'" = {e+p) u^'vT - pg^"" , (3) 

where s,p and u^^ are, respectively, the proper energy density, pressure and four- 
velocity of the fluid clement. Eq.(2) should be complemented by the continuity 
equations for conserved currents such as the baryon number, 

d'' (nuf,) = 0, (4) 

and also the equation of state which describes the thermodynamical properties of 
the matter. The equation of state may be specified in the form of 

£ = e (n, s) , (5) 

where n and s are the baryon number and the entropy densities, respectively. For 
ideal fluid case, we can also derive the conservation of entropy, 

d>' {su^) = 0. (6) 
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from Eqs.(2,3) and (4) using the Gibbs-Duhem relation, 



dp = jidn + Tds, 



(7) 



which can be uniquely determined once Eq.(5) is specified, together with the first 
law of thermodynamics, de = Tds + fidn and the Euler relation for extensivity 
e + p = Ts + /in. 

These hydrodynamical equations can also be derived from the action principle^. 
They can be solved numerically to give the space-time development of all the nec- 
essary thermodynamical variables and the fluid velocity^^'^. 



3. Some Open Problems of Relativistic Hydrodynamics 

As mentioned in Introduction, the ideal hydrodynamical description for the dy- 
namics of hot and dense matter observed in RHIC experiments works amazingly 
well, especially for the collective flow. The success of the approach indicates that 
the very early equilibration of the partonic gas takes place. The nature of the QGP 
seems rather a strongly interacting fluid (sQGP) than a ideal free parton gas^. On 
the other hand, there exist several open problems in the interpretation of data in 
terms of the hydro model. These questions require careful examination to extract 
quantitative and precise information on the properties of QGP. Furthermore, one 
should note that several very different physical scenarios within the hydrodynam- 
ical approach (e.g. the use of continuous emission, sudden freeze-out and rQMD 
final state cascade, or some drastically simplified hydro models such as blast wave 
solution, etc.) can give rise equally good results in reproducing the observables with 
suitable choices of parameters. In a way, one may say that the hydro signature is 
"robust" , but on the other hand, this could be a synonym of "insensitive" . In order 
to extract more precise information on the properties of the quark-gluon plasma 
and the mechanism of phase transition in the bulk QCD matter, we should clarify 
these points and refine the physical parameters used in these approaches. 

With respect to this point, one should remind a important point of the collective 
flow. When we use the equation of state to describe the pressure as function of 
density, the hydrodynamical equations are meaningful only if the approximation of 
local thermal equilibrium is valid. On the other hand, as far as Eq.(2) is concerned, it 
is nothing but the local conservation of energy and momentum density. Therefore, it 
may be possible that a system which is completely out of local thermal equilibrium 
can also manifest a flow pattern depending on the origin of the force. For example, 
as a well-known example, we may recall that the Schrodinger equation 



dtp „2 , T. , 
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can also be written just as if the form of hydrodynamics, 




_ + V • (nu) = 0, 




comes from the quantum uncertainty principle and serves as the pressure, but 
nothing to do with the thermal equilibrium. In this example, the "pressure gra- 
dient" comes from the quantum uncertainty principle. That is, a flow phenomena 
does not necessarily indicates the validity of the local thermodynamical equilibrium. 

In the following, we list several crucial points of the hydrodynamical approach 
to the relativistic heavy ion reactions. 

• Initial conditions 

When two nuclei collide, a huge number of partonic degrees of freedom are ex- 
cited. At the first instant of the collision, these partonic degrees of freedom are far 
from thermal equilibrium. To obtain the initial condition suitable for the hydro- 
dynamical scenario, we should understand how these partonic degrees of freedom 
reach the local equilibrium state, forming the quark-gluon plasma. It seems that 
the time scale to achieve the thermal equilibrium in partonic cascade calculation 
within reasonable values of parton cross section is much larger than the time scale 
required for the hydrodynamical models. Thus the success of hydrodynamics cast 
a very interesting question of how such an early thermalization can be attained 
in the initially created QCD partonic excitations. Several new concepts have been 
proposed. For example, as the mechanism of quick isotropization from the initial 
momentum distribution of partons (which is predominantly longitudinal), the so- 
called Weibel instabilities known in plasma physics is shown to be effective also in 



However, the isotropization itself is a necessary but not a sufficient condition for 
the thermal equilibrium. With respect to the mechanism of thermalization, another 
interesting notion, called pre-thermalization, is recently developed. In a scalar field 
theoretical model, it is shown that when the system is set to an excited state, the 
validity of equation of state (the functional relation between energy density and 
pressure) is attained well before the real thermal equilibrium is reached^^. With 
respect to this point, we have discussed a possible relation between the mechanism 
of early thermalization to the non-extended statistical mechanics^^. Furthermore, 
the approach from the color glass condensate may provide an initial momentum 



the case of QCD 



,11,12,13 
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distribution easy to achieve a quick thermalized state of partons because of satu- 
ration mechanism. Several investigations based on the CGC approach to calculate 
the initial energy distribution have been carried out^^. 

• Event-by-event fluctuations. 

Even the mechanism of thermalization is fast enough, the initial condition at- 
tained by high energy nuclear collision event is far from smooth as shown by the 
event generators such as NeXUS code^"'. Furthermore, this initial configuration fluc- 
tuates collision by collision very largely, even for central collision of heavy nuclei. 
It has been pointed out that the effect of even-by-event fluctuations due to the 
different initial conditions are crucial for quantitative studies of observables such as 
elliptic fl-Ow^®. Since the dynamical evolution of hydrodynamics is extremely non- 
linear, the event average of any hydrodynamical obserbable is quite different from 
the value calculated from the averaged smooth initial configuration, which is the 
hydrodynamcal calculations commonly used. The fluctuation of V2 coefficient must 
contain the dynamical information of the state formed in the early stage of the 
collision. 

• Finite size effect 

Strictly speaking, the hydrodynamics is essentially the zero mean-free path (or 
correlation length) approximation for the microscopic degrees of freedom compared 
to the system size. For the nuclear collisions, this is not a trivial approximation. 
As we know that the nuclear ground state has finite surface thickness (of the order 
of 3 fm) but if we try to describe the density distribution in terms of zero mean 
free-path hydrodynamics, we would have zero surface thickness for the ground state 
mass distribution. The parameter which characterizes the degree of approximation 
with respect to the finite size of the system would be 

^micro 

1^ 

-^Hydro 

where Xmicro is the typical mean-free path (or the correlation length) of the con- 
stituent particles and Xnydro is the hydrodynamical inhomogeneity scale, which is 
typically 



A 



Hydro 



p 



(12) 



Hydrodynamics is valid for a; ^ 1, so that when the local inhomogeneity scale in 
the system is of the order of the mean-free path, the hydrodynamics should break 
down. 

In order to take into account for the finite x values, as relativistic extension of 
the Weizsacker term, we may introduce the gradient term in the action of hydro- 
dynamics as 



-J 



d^x 
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where a is a quantity related to the surface tension in the static hmit and is 
the four-gradient projected to the hyperplanc perpendicular to the velocity, 

where P^" = (5^" — u'^u") is the projection operator to the hypersurface perpen- 
dicular to the four- velocity u^^. The general energy-momentum conservation can 
be obtained from the Nother theorem and the resulting equation of motion (for 
constant a case) becomes 

d^f^"" = -d^ ((7(V'^n)(5^n) - anu'^V^n - ^g^''iV^n){drn)^ 

where 

f^"" = suV - P^^p 

with 

p = p + andi,{y^n). 

Unfortunately, the above equation becomes dynamically much more complicated 
than the usual hydrodynamics. In any case, the finite size effect on the local energy 
density is an important question, in particular with respect to the properties of 
phase transition. A further investigation on this direction is now in progress^^. 

• Viscosity effects. 

The comparison of flow phenomena to the ideal fluid calculation indicates that 
the emergence of new state of matter which flows like an ideal-fluid at RHIC ener- 
gies. This idea that the QGP behaves a real ideal fluid raised an interesting question 
when the viscosity for the strong coupling limit of 4D conformal theory obtained 
from the supcrsymmetric string theory found to be a very small^^. On the other 
hand, Hirano and Gyulassy argue that this is due to the entropy density of the 
QGP which is much larger than the hadronic phase^. To be precise, a quantitative 
and consistent analysis of the viscosity within the framework of relativistic hydro- 
dynamics has not yet been done completely. This is because the introduction of 
dissipative phenomena in relativistic hydrodynamics casts difHcult problems, both 
conceptual and technical. Several works have been done in this direction'^^. We will 
discuss later this point more in detail in the next section. 

• Final hadron spectra from the hydrodynamical model. 

To analyze the physical observables in terms of the hydrodynamical scenario, 
we have to construct the particle spectra from the hydro solution. As the hydro- 
dynamical expansion proceeds, the fluid becomes cooled down and rarefied, thus 
leading to the decoupling of the constituent particles. At this stage, these particles 
do not interact any more and free-stream to the detector. Long-lived resonances 
and other unstable particles may decay on the way to the detector after this instant 
of decoupling phase. In the standard hydrodynamical models, one introduces the 
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concept of freeze-out, which assumes that particle emission occurs on a sharp three- 
dimensional surface (defined for example by the local temperature, T{x,y,z,t) = 
constant). Before crossing it, particles have a hydrodynamical behavior, and af- 
ter they free-stream toward the detectors, keeping memory of the conditions (flow, 
temperature) of where and when they crossed the three dimensional surface. The 
Cooper-Prye formula '^^ gives the invariant momentum distribution in this case 



da^ is the surface element 4- vector of the freeze out surface a and / the thermal 
distribution function of the type of particles considered. The space-time depen- 
dence comes from those of thermodynamical parameters, such as temperature and 
chemical potential. This is the formula implicitly used in all standard thermal and 
hydrodynamical model calculations. 

Although simple and elegant, this sudden freeze-out is not only an idealization 
but also contains some problems such as conservation of energy and momentum, 
negative flux and artiflcial entropy production*. To remedy this approach several 
approaches have been proposed such as the use of URQMD code coupled to the 
final state of the hydrodynamics^^'^^. However, usually these calculations are rather 
complex and time consuming in practice. 

The most important effect which should effect the form of particle spectra is 
that not every particles are emitted from the same hypersurface specified by a 
unique temperature. The continuous emission approach^^ still uses the equilibrium 
momentum distribution of particles but they can be emitted continuously during the 
hydrodynamical evolution of the system. It also takes into account the absorption 
effects while the emitted particle from the inside traverses the surrounding hadronic 
matter. This accounts for the finite size of the system in the final particle spectra. 
In contrast to the usual sudden freeze-out, it is found that the final hadrons can be 
emitted from a broad range of temperatures. It should be emphasized that the two 
extremely sharp temperature surface of sudden freeze-out and almost fiat 

distribution of temperature of the continuous emission scenario give equally good 
description of the observed spectra and fiow^. 

4. Viscosity and Causality 

To extract more precise quantitative conclusion on the ideal nature of the QGP 
fluid, we should study the effect of dissipative processes on the collective flow vari- 
ables. However, a covariant theory of dissipative phenomena is not trivial at all. 
We know that the diffusion equation is not covariant, and the equation is parabolic 
so that the propagation of signal has an infinite velocity. Landau introduced the 
dissipative effects in the relativistic hydrodynamics, but it is known that the formal- 
ism of Landau^ of relativistic viscous fiuid still leads to the problem of acausalily. 
Relativistic covariance is not the sufficient condition for a consistent relativistic 
dissipative dynamics. To cure this problem, the second order thermodynamics was 




(18) 
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developed by Israel, Stewart and Miller^^. However, this theory is too general and 
complex, containing many unknown parameters which make difficult the applica- 
tion of the theory to practical problems. Furthermore, the theory contains the third 

order time derivatives, introducing additional difficulties of initial value problem 
and numerical procedure. In our opinion, this is not the unique approach to the 
consistent relativistic dissipative hydrodynamics. Here, we show that an alternative 
theory which satisfies the minimum conditions of covariance and causal propagation 
of signal, and in the small viscosity limit, it recovers the usual ideal hydrodynamics. 

To this, we first note that the problem of acausal propagation in usual diffusion 
equation, 

can be cured by the introduction of the relaxation time as 

d'^n dn 

Trelax-Q^ + ^1, (20) 

converting the parabolic equation to the hyperbolic equation. For a suitable choice 
of parameters, Treiaxi we can recover the causal propagation of the diffusion flux??. 
Physically this can be understood as following. In general, the diffusion equation is 
the combination of two equations, 

• Continuity equation, 

^ + V.J=G, (21) 

• Generation of the irreversible current due to the thermodynamical force, 

and within the linear response of the system, we have 

3 = -CVn 

where F is the thermodynamical potential and C, is in general a function 
of thermodynamical quantities, but here we assume to be constant for the 
sake of simplicity. 

The origin of acausality resides in Eq.(22) than the continuity equation, (21), 
since this means that the action of the force V {5F/ 5n) immediately generates the 
physical current j. We may think of relaxation phenomena for this process. This 
can be done phenomenologically by introducing the retardation function 

G {t, t') = ^^e-(*-*')/^-'- (24) 

and rewrite Eq.(22) as 

ft 

G{t,t')\/F {f,t')dt'. 
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In the limit of Treiax — 0, we have G {t, t') 5(t — t') so that this recovers the 
original equation, (22). Now 

J = VF j. 

OZ '^relax '^relax 

Taking the time derivative of Eq.(21) and substituting the above equation, 
d'^n „ / 1 „^ 1 



= V- Vf 



'^relax ^relax 

1 ^ V^F 



'^relax \ 



obtaining Eq.(20). Eq.(20) is usually referred to as telegraphic equation. For more 
microscopic derivation of telegraphic equation, see 28,29,30 

Wc use the above mechanism to derive the relativistic dissipative hydrodynamics 
to transform the Landau formulation to satisfy causality. For this purpose, first let 
us review briefly the essential part of the Landau derivation of the relativistic 
dissipative hydrodynamics. Landau requires the conservation laws. 



d^T^^" = 0, (27) 
S^AT^ = 0. (28) 

In the presence of dissipative phenomena, the energy momentum tensor T'^" and 
the baryonic current N'^ are not given by Eqs.(3) (4) anymore, but instead, 

T^"" = eu^u" - P"" (p + n) + TT^'' , (29) 

N" = nu" + 1^", (30) 

where 11 and tt^'' are the bulk and shear viscous stresses, respectively and is the 
diffusion current of baryon number. For these, we require the constraints, u^tt'^" = 
and UnV^ = 0. With these terms, of course the conservation of entropy Eq.(6) is 
not valid anymore, instead we have 

(sw" - au") = ^ i-P'^-'n + tt"^) d^u^ - u^d^a, (31) 

where a = n/T. Landau identifies the term 

0-^ = su^' - ai/^ (32) 

as the entropy current and requires the positive definiteness of the right-hand side 
of Eq.(31) as 

L^-pf^-U + nnd^u,,-iy''d^a>0, (33) 
obtaining the following expressions for the viscous terms, 

n = Cdau", TT^v = r]dnU„, = -Kdf^a (34) 
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and introduce the projection operator to satisfy the orthogonality condition for the 
viscous term to the fluid velocity as 

where pf^°"^>^ is the double symmetric tracclcss projection, 

p^iauj} _ \ ^p^a pvfj _|_ pua pL'v'j ^ pnu pajS (36) 

2 3 

As mentioned, this Landau scheme leads to the acausal propagation of thermal 
current. However, at this moment, generalization of these equation in order to 
obtain to hyperbolic equations is self-evident. We introduce the retardation integral 
in each viscous term before the projection, 

n = - dr'G (r, r') C9„«« (r') , (37) 
#'''^ = J\T'G{T,T')d''u'' {t'), (38) 

i"" = - j dr'G (r, r') kQ^Q! (t') , (39) 

where r is the local proper time. These integrals are equivalent to the differential 
equations 

n = -ca„«« + 7^ (40) 

jr^"" = rjOf^u" - 7-— (41) 
ar 

= -Kd^a + 7^ (42) 
ar 



with 



is the total derivative with respect to the proper time. These equations, after the 
projection (Eq.35), can be compared to the Israel-Stewart form, 

n = -C (^daU- - /^o ^ - "oSa'^") (44) 
daU0 - - a,daU0\ (45) 

= -npi'- \d^a-(3^—+ a.djl + j . (46) 

They are similar, but we can see that the Israel- Stewart form contains the more 
general linear combination of the second order variables. However, the most im- 
portant and essential difference is that, in our formalism, the projection operators 
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P^" enter after the integration of these equations. All the dissipative terms are 
expressed explicitly in terms of independent variables of the usual hydrodynam- 
ics. In the Israel-Stewart form, due to the projection operators, integral form as 
Eqs. (37,38,39) can not be obtained explicitly. 

Our approach has several advantages to the Israel-Stewart formalism. First of 
all, we keep the simple physical structure of Landau formalism, but at the same 
time we cured its causality problem by introduction of the relaxation integral, 
Eqs. (37,38,39). Use of these integral expressions also eliminates the problem of 
higher derivatives in time. We only need the past values of the independent vari- 
ables to solve numerically. This also eliminates the problem of extra initial con- 
ditions, too. The integral expressions are easy to be evahiated when we use the 
Lagrangian coordinate system such as SPHERIO^"'^. Because of the simple form 
of viscous terms, incorporation of these equations to the realistic hydro-code such 
as SPHERIO is relatively easy. A work on this line is in progress. 

5. Summary 

The hydrodynamics is found to be a very successful tool for the description of 
the relativistic heavy ion collisions. However, from the quantitative point of view, 
the present hydrodynamical approach still contains many uncertainties and also 

some conceptual problems. In this paper, we call attention to these questions, and 
discussed some of problems in detail. In particular, we propose an alternative the- 
ory to the Israel-Stewart second order thermodynamics, where viscous terms are 
given by the integral expressions which take into account of the relaxation time. 
In this way, the problem of causality is avoided and at the same time a simple 
physical structure of Landau formulation has been kept. Other questions as early 
thermalization, finite size effects, fluctuations in initial conditions, etc. should also 
be studied more in detail When these questions are clarified, we will have much 
more detailed knowledge about the dynamics and properties of the new states of 
strongly interacting matter. It is also quite important to look for signals of genuine 
hydrodynamics, such as shock wave propagation^^. 

This work is dedicated to Prof. W. Greiner, one of the pioneers of the sub- 
ject, on the occasion of his 70th birthday. The authors express their thanks to 
H. M. Kiriyama, F. Grassi, Y. Hama, C.E. Aguiar and E. Fraga for stimulating 
discussions and kind help. T. Kodama is also grateful for the kind hospitality of 
Prof. H. Stocker during his stay in Frankfurt. This work was partially supported 
by FAPERJ, FAPESP, CNPQ and CAPES. 
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